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The recent passage of Jupiter by the quasar QSO J0842+1835 at a separation of 3.7 arcminutes 
on September 8, 2002, combined with recent advances in interferometric radio timing, has allowed for 
the first measurement of higher-order post-Newtonian terms in the Shapiro time delay which depend 
linearly on the velocity of the gravitating body. Claims have been made that these measurements 
also allow for the measurement of the propagation speed of the gravitational force. This conclusion 
disagrees with recent calculations done in the parameterized post-Newtonian (PPN) model, which 
find no dependence of the velocity-dependent terms in the time delay on the "speed of gravity" to 
the stated order. Here, to test out these claims and counterclaims, we calculate the time delay in 
the limit of an instantaneous gravitational force, and find that the velocity-dependent terms are in 
complete agreement with previous PPN calculations, with no dependence on the speed of gravity. 
' We conclude that the speed of gravity cannot be determined by measuring these terms in the Shapiro 

, time delay, and suggest a reason why other groups mistakenly came to the opposite conclusion. 

I. HISTORICAL BACKGROUND 

Almost 40 years ago, I.I. Shapiro pointed out that the time delay which results from light appearing to slow down 
as it passes through a gravitational potential could be measured within our solar system [1]. The measurement 
was done three years later, and provided a new test of the theory of general relativity (GR) [2]. Recent advances 
\ have now allowed radio astronomers to measure higher order post-Newtonian (PN) terms in the Shapiro time delay, 
using Very-Long Baseline Interferometry (VLBI). A recent passage of Jupiter by the position of a bright quasar at 
a separation of only 3.7 arcminutes has allowed for some of these high-order terms to be measured by Fomalont and 
Kopeikin, specifically those that have to do with the transverse motion of Jupiter perpendicular to the line of sight 
^C) \ to the quasar [3]. In a series of papers, Kopeikin has argued that the terms in the Shapiro time delay which depend 
upon the velocity of the gravitating body also depend upon the propagation speed of gravity [4,5]. Specifically, in 
"*^h | Ref . [4] , he derives an equation for the time delay is the form 
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where td is the time delay, m a and v a the mass and velocity of the a'th gravitating body, respectively, c g is the 
^ , propagation speed of gravity, f a the separation vector between the gravitating body and the observer, r a its magnitude, 
C is an integration constant, and the vector K (called N by Kopeikin) is given by 
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where k is a unit vector pointing in the direction of the photon's path. In what follows, for clarity, we will refer to 
equations in Ref. [4] with the author's initial, thus our Eqs. 1 and 2 correspond exactly to Eqs. (K22) and (K23). 

Recently, Asada [6] has taken issue with the interpretation of the velocity-dependent terms in the time delay 
for models which assume that c g — c. Will goes further, evaluating the time delay for models in which c g c, and 
concludes that there is an error in Kopeikin's derivation of the time delay [7] . Using the parameterized post-Newtonian 
(PPN) framework (see Ref. [8] for more details), he concludes that the proper time delay is given by 

*<* = - 2 E f 1 - — ) ln (^ - R • ? «) + C ' ( 3 ) 



where 



K = k k x (v x k). (4) 
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Using his initials for clarity as well, our Eqs. 3 and 4 correspond exactly with Eqs. (W35) and (W36), where we have 
set several of the PPN parameters which do not affect the final result equal to the values predicted by GR. 

In what follows, we will repeat the calculations of Kopeikin and Will, but in a model where the gravitational force 
is instantaneous, i.e. c g — > oo. We identify where the two methods agree under this assumption, and where they 
find different mathematical expressions. In short, we find that Will's result is correct given his assumptions regarding 
the PPN framework. Kopcikin's result is inconsistent with his assumptions, revealing that his derivation contains a 
mathematical error. We identify the likely cause of it in the course of our derivation. 



II. CALCULATING THE TIME DELAY 



For the following calculations, we make use of the following assumptions. A photon travels along the x-axis of our 
coordinate system. It passes a body of mass m a moving with velocity v a = [v\\ , v±, 0) along a straight line. The origin 
of the time coordinate is defined such that the instantaneous separation vector between the gravitating body and the 
photon is perpendicular to the photon's path at t = 0. This does not correspond exactly to the moment of closest 
approach between the photon and the gravitating body, which can easily be calculated to be first-order in v a /c. Using 
these velocities, we find that the position of the of the gravitating body is given by x a (t) — (vnt, v±t + yo, zq) and that 
of the photon by x p (t) = (ct, 0, 0). Our observer is placed far from the gravitating body at some distance x along the 
photon's path, at position x = (x o ,0,0). As these calculations are traditionally done by placing an observer at the 
barycenter of the system, we should technically include at least one other gravitating body, but we see immediately 
that it will have no effect to lowest order on the measured time delay. 

To calculate the time delay for the photon passage, we note that since we can make the gravitating mass sufficiently 
small, we can ignore all effects associated with the deflected path of the photon, which are of higher order in v a /c. 
The equations of motion for the photon's path as a function of its wave vector k^ = dx^/dt = (1, k l ) are completely 
determined from a description of the spacetime metric, and the constraint equations for a null geodesic. These 
constraint equations are given by both Eq. (Wll) and Eq. (K14), in complete agreement, in the form 

<h l ' lp + k^(r^ - fcT°J = (5) 
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ds 

— J = goo + 2g Qi k l + gijW = 0. (6) 

While Kopeikin and Will agree on the geodesic equation for the photon, they disagree on the form of the metric. 
Will, using the standard parameterized post-Newtonian (PPN) formalism, assumes that the metric, Eq. (W9), takes 
the form 

g 00 = -l + 2U (7) 
90i = -4U, (8) 
g ij = (l + 2U)6 ij , (9) 

where we have set the PPN parameters 7 and ai equal to the GR values 7 = 1, ct\ — 0, and the retarded Lienard- 
Wiechert mass and momentum potentials are given by Eq. (W10) as 



U(t,x) = ^^ (10) 
c 2 r e {t,x) 



where we define r e to be the "effective" distance of the retarded potential, such that r e (t,x) = \x — x a (t')\, where the 
"emission" time of the gravitational force must satisfy t — t' = r e /c g , where c g is the speed of gravity. 

Kopeikin's metric is slightly different, in ways that deserve some clarification. He introduces a parameter r, which is 
used to describe all quantities involving gravitation in his framework. It is defined such that c g r = ct. Unfortunately, 
the gravitational equations treat this parameter as a physical time for gravitation. Thus, the wave equation for linear 
metric perturbations, Eq. (K5), reads 

^^2+ V 2 ) r»(T,x) = -l^e-(r,x), (12) 
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where 7^" is the standard trace-reversed linear metric perturbation, and M1/ is his modified stress-energy tensor. We 
note that for this equation to give a traveling wave solution moving at speed c g , we must interpret r as the physical 
time, not t. We believe it is confusion between these two quantities which leads to errors in his final conclusions. In 
any case, Kopeikin modifies not only the stress-energy tensor but the metric perturbation as well, yielding a metric 
that agrees with Will's (if we view r as the physical time-variable) but for a different momentum potential, given by 
Eq. (K10) as 

T/ ,. Gm a (v a )i 

m > x)= JWtfj- (13) 

In the calculation that follows we will see that the differences in the two metric formulations lead to different 
values for the gravitomagnetic drag experienced by a photon passing by a moving body. We will also see that Will's 
calculation is completely consistent, whereas Kopeikin's final answer is inconsistent with the photon propagation 
equation. 

We start our calculation by deriving the magnitude of the photon's wavevector in our metric, Eq. 8, noting that 
Kopeikin and Will differ on the value of the momentum potential V,. We denote the magnitudes of the wave- vector 
k l = (k, 0, 0) and the momentum potential Vi by k and V, respectively. The x-component of the momentum potential 
is denoted Vj| . From Eqs. 6-9, we find 

dr AV W + \l l - 4JJ2 + UV \\ 
(-1 + 2U) - 8V {l k + (1 + 2U)k 2 = - k ee = V__ L. (14) 

Working only to 1.5PN order, we throw out all terms involving U 2 or V 2 , and find 

dr 1 + 4VIi 

k =ir = tw~ (1 - 2C,+41 * <»' 

To find the time delay, we integrate this expression over the path length, finding for the time delay 

«-</*>-*-[/ ^T^nd " * ~ / dJ ? {W ~ iV ' ] ' (16) 

where to is the time required to travel a distance Sx in the absence of gravitating bodies. The two terms here have 
entirely different meanings, but both are easily understood. The first, oc / Udx p , which appears for static gravitational 
sources as well, is the simple geometric time delay which results from photon's traveling through any gravitational 
potential. The second term, cx J V\\dx p is the gravitomagnetic contribution to the time delay. In simplest terms, it 
can be thought of as a "frame-dragging" term, whereby the moving object pulls the photon along it's path. At this 
point, it is useful to compare our result with those of Will and Kopeikin. To do so, we make use of the ratio between 
the momentum potential and the mass potential, which is constant throughout space and time in both formalisms so 
long as the gravitating body moves at constant velocity. From Eqs. 10 and 11, we see that in Will's formalism that 
Vi = [/i^ali, whereas from Eqs. 10 and 13, Kopeikin's method yields Vi = JJ^^i. We find, respectively, 

C Cg 

t d ~2(l- J (Will) (17) 

~ 2 ^1 - J (Kopeikin), (18) 

which correspond with Eq. (W16), and Eq. (K21). 

Calculating the time-delay integral is straightforward, so long as we calculate the effective distance correctly for 
our chosen value of c g . To simplify matters, and emphasize the difference between the frameworks used by Will and 
Kopeikin, we will evaluate the time delay for a model with instantaneous gravitational propagation, i.e., c g — > 00. In 
Kopeikin's notation, this represents the limit e — > 0, which he refers to as the Newtonian limit. Given the trajectories 
of the photon and gravitating body stated above, we find that the effective radius is given as a function of time as 

r e (t)= [( C -«||) 2 i 2 + Kt + yo) 2 + ^]°' 5 

= [(c 2 - 2cw|| + v 2 a )t 2 + 2y v ± t + y 2 + z 2 ] °' 5 . 
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To simplify the expression, we note that at the moment of closest approach, at t m in = —yov±t/[(c 2 — 2cv» + v^)t 2 }, 
the distance reaches its minimum value of T m £ n = 

r e (t m in)- At all other times, the effective distance is given by 
r e (t) = fL + IWi-t mm ))P, (19) 

where u = k — v a /c, k is the unit 3- vector pointing along the photon's path, and we know that f min -u = 0. Note that 
the appearance of c in these expressions has nothing to do with assumptions about the speed of gravity, it merely 
represents the speed of the photon when we calculate relative distances and velocities. 

The calculation of the time delay will allow for a direct test of Will and Kopcikin's results, since both derive the same 
integral for the time delay, albeit with different prefactors, as we found in Eqs. 17 and 18. The time delay integral is 
analytic, since we can integrate along the unperturbed path of the photon. Changing variables to dt — dx p /c, shifting 
the time axis to a new time variable offset by t min , and defining u = \u\, we find 

f dt „ f dt „ f dt 
td ~ / 7777T ~ Gm - / ~m ~ Gm - / / 2 1 2 m (20) 

J U{t) J r e (t) J V r mm + U C t 



~ Gm a (- j H^r 2 min + uW f - uctf], (21) 

where tf is the arrival time of the photon. Note that in deriving this expression, we made no assumptions about 
the geometry of the system at the moment of closest approach, which Kopcikin discusses in his appendix. In what 
follows, we will only use the fact that vectorial velocity difference between the photon and the gravitating body is 
perpendicular to the separation vector at closest approach. 

We start with an analysis of the prefactor 1/uc. Calculating this to first order in w a = v a /c shows us that 

■1 = ^— = - 1 = - 1 (22) 

uc c\k-w a \ r Jk.k_ok., 1 ,_+, lT _., 1 t_ c^l-2w ll+ w 2 



cy k ■ k — 2k ■ w a + w a ■ w a 

- n 1 , ~ -(1 + w\\) = -(1 + v\\/c) = -(1 + ^). (23) 
c(l —w\\)c c 11 c c 

This multiplicative factor is derived directly from the integral of the gravitational potential along the line of sight, and 
has absolutely nothing to do with the speed of gravity, which we are assuming here to be instantaneous. It is essentially 
a Dopplcr correction, representing the increased (decreased) time that the photon spends in the deepest part of the 
gravitational potential because of the gravitating body's parallel (anti-parallel) velocity component. Comparing to 
Eq. 17, we see that the "Doppler" term is half the magnitude and opposite in sign from the gravitomagnetic correction 
term. On the other hand, this result is at odds with the gravitomagnetic term derived by Kopcikin in Eq. 18, which 
contains c g in the denominator. We will discuss this error below. 

Moving on to the logarithmic term in the integral, we define the separation vector r(t) = x p (t) — x a {t) = r min + ctu 
whose length is r e (t) — \/rf nin + u 2 c 2 t 2 . Additionally, we note that 

k ■ r(t) = k ■ r min + ct(k ■ k - k ■ w a ) = k ■ r min + ct(l - ^_^). (24) 

c 

This expression is similar to the term uct which appears in logarithmic term of the time delay, but not equal to it at 
the required order because of the appearance of a term k ■ f min . However, we see that 

T'min ^ ► fmin ' k fmin ' 0- (25) 

Denoting by r\\ and r± the components of f m in parallel to «|| and v±, respectively, we find, to lowest order, 

r±w± = r=|| (1 - W|| ) ~ fj| = f min ■ k. (26) 

Thus, we can subtract off the term f m in ■ k in Eq. 24 by taking the product of the transverse velocity of the gravitating 
body with the minimum separation vector. It is trivial to deduce that k x (v a x k) = v±. We see though, that to first 
order in v a /c, 

W± ■ f(t) = W± ■ fmin + ct(w± ■ k - Wj_ ■ w a ) (27) 
~w ± r ± , (28) 
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where the second term vanishes from orthogonality, and the last because it is second-order in v a /c. Combining these 
results, we conclude 



Plugging Eqs. 21, 23, and 30 into Eq. 17 we recover Eqs. (W35) and (W36). Plugging the same three equations into 
Eq. 18 does not recover Kopeikin's Eqs. (K22) and (K23), which have the speed of gravity c g where we have the speed 
of light c. We conclude that his calculation is erroneous, since the time-delay integral we calculated had nothing to 
do with the speed of gravity, which we took to be instantaneous. It is also independent of all metric quantities other 
than the gravitational potential which appears in goo . Why then does he find different numbers than ours? 

A clue is provided by the fact that our results (as well as those of Will, whose results ours reproduce) can be brought 
into agreement with Kopeikin's if we replace v a , the velocity of the gravitating body, with v a c/c g . In his notation, 
this is the transformation v a — > ew a . We believe the error he makes is to use r = et as the physical time for the 
gravitating object, declaring that its physical velocity is given by dx a /dr — v a , while keeping t as the time variable 
used to describe the photon's motion, as in Eq. (K13). If we convert the velocity he uses for the gravitating object 
from r-based coordinates to t-based coordinates, we find immediately that v a {t) = dx a /dt = dx a / dr{dr / dt) = ev ai 
which corresponds to the unphysical and erroneous velocity found in his calculations. 



We have repeated the calculations performed by Will and Kopcikin, in a model where the speed of gravity is 
instantaneous. We find that first-order terms in the time delay resulting from the motions of gravitating bodies are 
independent of the speed of gravity, and follow the form written down by Will in every respect. We find in addition 
that Kopeikin's equation for the time delay, which asymptotically approaches the time-delay for static bodies in the 
limit of an infinite speed of gravity, is wrong. We conclude that since the first-order velocity terms in the time delay 
are independent of the speed of gravity, Fomalont and Kopeikin's high-precision measurements of the time delay of 
light from quasar QSO J0842+1835 passing by the edge of Jupiter on September 8, 2002 [3], while a truly impressive 
observational feat in radio astronomy, provided a measure of the speed of light only, not the speed of gravity as was 
claimed. 

The author wishes to thank P. Grandclement and E. Bertschinger for helpful discussions. He also wishes to thank 
S. Kopeikin for clarifying discussions about the details of his results. 
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